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We calculate the thermal photon transverse momentum spectra and elliptic fiow in y^sjvjv = 
200 GeV Au+Au collisions at RHIC and in ^/s^ = 2.76 TeV Pb+Pb collisions at the LHC, using 
an ideal-hydrodynamical framework which is constrained by the measured hadron spectra at RHIC 
and LHC. The sensitivity of the results to the QCD-matter equation of state and to the photon 
emission rates is studied, and the photon V2 is discussed in the light of the photonic pr spectrum 
measured by the PHENIX Collaboration. In particular, we make a prediction for the thermal photon 
Pt spectra and elliptic flow for the current LHC Pb+Pb collisions. 



PACS numbers: 12.38.Mh,25.75.Cj,25.75.Ld 

I. INTRODUCTION 

Experimental data at the Relativistic Heavy Ion Col- 
lider (RHIC) and now also at the Large Hadron Collider 
(LHC) have shown compelling evidence of strongly in- 
teracting QCD-medium production in ultra-relativistic 
heavy-ion collisions. The measured transverse energies, 
transverse momentum (pt) spectra and, in particular, 
the significant azimuthal anisotropy (elliptic flow) of 
final-state hadrons suggest together that partonic QCD 
matter, quark-gluon plasma (QGP), is formed in these 
collisions. 

In the collision of two nuclei, the azimuthally 
anisotropic overlap region sets preferred directions for 
the transverse flow. In hydrodynamical models, pres- 
sure gradients turn the spatial anisotropy of the produced 
hot matter into a flow anisotropy, which is transmitted 
into the momentum distributions of measurable final- 
state hadrons at the decoupling of the system. However, 
the hadronic measurement reflects the flow (and temper- 
ature) conditions only in the freeze-out region where the 
hadronic interactions cease. 

In comparison with partons, photons interact only very 
weakly with the QCD matter and thus a photon emitted 
from the medium most likely escapes from the system 
without interacting. This is seen also in the measure- 
ments where photons do not show a similar suppression 
as hadrons when we move from proton+proton (p+p) 
to nucleus-f nucleus {A A) collisions [l| . Since photons 
can escape from the medium without interacting, they 
carry information about the system at the time of their 
production. 

Ultra-relativistic heavy-ion collisions are particularly 
interesting in regard with direct photon production, since 
relative to p-l-p collisions there are different types of nu- 
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clear effects at work as well as a number of important 
further sources for photons. In p-l-p (and also in \i-\-A) 
collisions the direct photons are prompt photons orig- 
inating from the primary hard interactions of partons, 
and fragmentation photons emitted by the primarily pro- 
duced high-pT partons [5]. In heavy- ion collisions (and 
also in p-|-A), both the prompt and fragmentation pho- 
tons at high-pT are subjected to nuclear effects in the 
parton distribution functions of the colliding nuclei (see 
e.g. [1, 01 for the quantification of these effects and their 
uncertainties). The fragmentation photon component is, 
however, expected to be suppressed due to the quench- 
ing of partonic jets in QCD matter m A + A collisions. 
In addition to this, m A + A collisions the jet(parton)- 
matter interactions, i.e. the jet-photon conversion |g| 
and collision-induced photon emission from high-energy 
partons can produce photons which are important in the 
mid-pT region [Z^j)] . Finally, the hot medium itself emits 
thermal photons, which are expected to be important in 
the few-GeV region and below, as discussed in the hy- 
drodynamical studies of Refs. i7:, 8, 10-13]. 

In heavy-ion collisions, it is very difficult to distinguish 
between the different direct photon sources. In addition, 
there is a huge decay-photon background to deal with. 
The elliptic flow of direct photons could, however, shed 
more light on the interplay of the various photon pro- 
duction sources which differ from each other as follows: 
At high-pT (above ~ 5 GeV at RHIC), where prompt 
photons dominate Q , and where the fragmentation pho- 
tons are more suppressed in the out-of-plane direction 
(perpendicular to the impact parameter), the photonic 
V2 should be positive but very small Q. The jet-medium 
interactions in turn increase the photon production most 
strongly in the in-plane direction (parallel to the impact 
parameter), thus causing a negative V2 contribution at 
mid-pT [s|- The thermal photon production is affected 
by the hydrodynamical transverse flow itself, so that pho- 
tons in the in-plane direction get a stronger boost. As 
shown earlier in Refs. (14i - il6j . this results in a positive 
elliptic flow for the thermal photons. Since the net contri- 
bution from other sources to photon V2 is expected to be 
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very small or even negative ^ , a large (hadron-like) pho- 
ton V2 measured in the few-pr region and below, should 
thus serve as a signature of thermal photon dominance. 
Since QCD matter is emitting photons throughout its en- 
tire evolution, measuring thermal photon px spectra and 
V2 would thus give important further constraints for the 
dynamics and properties of QCD matter. 

In this work we focus on computing the thermal photon 
Pt spectra and elliptic flow in ^/snn = 200 GeV Au-|-Au 
coUisions at RHIC and in ^/sWn = 2.76 TeV Pb-hPb col- 
lisions at the LHC, using an ideal- hydrodynamical frame- 
work which is constrained by the measured hadron spec- 
tra at RHIC and LHC. We study the sensitivity of the 
results to the QCD-matter equation of state (EoS) and 
to the photon emission rates. We discuss the photon V2 
in the light of the photoi iic p t spectrum measured by the 
PHENIX Collaboration |17H19| . In particular, we make a 
prediction for the thermal photon px spectra and elliptic 
flow for the current LHC heavy-ion collisions. Previous 
predictions for the thermal photon pT spectra in Pb-|-Pb 
collisions at the planned maximum cms-energy 5.5 TeV 
of the LHC can be found in [13 [2l|. 

II. THEORETICAL FRAMEWORK 
A. Centrality classes 



shown in [23|, we can get a good description of the 
pion spectra and the elliptic flow with this pQCD + sat- 
uration + hydrodynamics approach. For Au+Au colli- 
sions at y^SNN = 200 GeV the model gives an initial 
time To = 0.17 fni. For the ^s]^ =2.76 TeV Pb+Pb 
collisions at the LHC, we fix the initial entropy so that 
we reproduce the measured multiplicity ^25^. The initial 
time To = 0.12 fm is based on the EKRT- motivated fit 
done in Ref. 

To fix the initial transverse density profile in y^STvw = 
200 GeV Au-|-Au collisions, we do the following: In 
Fig. [T]we show, from Ref. [23|, the measured charged- 
particle multiplicity (divided by the number of partici- 
pant pairs) as a function of the number of participants 
calculated from the optical Glauber model^. Choosing 
the initial transverse density according to the binary- 
collision-scaled energy or entropy density (eBC, sBC), 
or wounded-nucleon-scaled energy or entropy density 
(eWN, sWN) as introduced in Ref. [l^, we compute the 
charged-particle multiplicity in the centrality classes ob- 
tained above. The initial entropy at 6 = in these four 
cases is kept fixed. We see that the sWN profile fits 
the measured centrality dependence quite well. We will 
therefore choose the sWN profile at RHIC and, for sim- 
plicity, use the same profile also for the LHC Pb-|-Pb 
collisions. 



Centrality classes for A + A collisions studied here are 
calculated using the optical Glauber model. For nuclear 
densities we use spherically symmetric Woods- Saxon pro- 
files with the the thickness parameter d = 0.54 fm and 
radii Rau = 6.37 fm and Rph = 6.49 fm. The total cross 
section ioT A + A collisions is calculated from 

= / d'b^ - J d^b(l-e-^-(^)'^W"„), (1) 

where Taa is the standard nuclear overlap function and 
crj^jy is the inelastic nucleon-nucleon cross section. We 
take crj^jv = 42(64) mb for ^/JWn = 200 (2760) GeV. 

The centrality classes are defined with impact param- 
eter ranges [6i,6i+i] so that for the centrality class of Ci 
we have 




The average impact parameter for each centrality class 
is calculated using the distribution da/(Pb as a weight. 
The average number of participants is calculated simi- 
larly. The obtained centrality classes, impact parameter 
ranges, average impact parameters and number of par- 
ticipants are listed in Table H] 

B. Initial states for hydrodynamical evolution 



C. Hydrodynamics and freeze-out 

To describe the spacetime evolution of the produced 
QCD matter, we solve the ideal-hydrodynamic equations 

= 0, (3) 





centrality % b 


range [fm] 


(fe) [fm] 






0-5 


0.00-3.35 


2.24 


346 




5-10 


3.35-4.74 


4.08 


289 




10-15 


4.74-5.81 


5.30 


242 




15-20 


5.81-6.71 


6.27 


202 




20-30 


6.71-8.21 


7.49 


153 


RHIC 


30-40 


8.21-9.49 


8.87 


102 




40-50 


9.49-10.6 


10.1 


64.4 




50-60 


10.6-11.6 


11.1 


37.5 




0-20 


0.00-6.71 


4.47 


267 




20-40 


6.71-9.49 


8.18 


128 




0-5 


0.00-3.53 


2.35 


375 


LHC 


0-20 


0.00-7.05 


4.70 


294 




20-40 


7.05-9.98 


8.60 


141 



TABLE I: Various centrality classes for Au+Au collisions at 
^^SNN = 200 GeV and for Pb+Pb collisions at y^sjviv ~ 
2.76 TeV, obtained via the optical Glauber model. 



For RHIC we use the EKRT saturation model ^ 
to fix the initial entropy in most central collisions. As 



^ Note that usually the number of participants quoted by the ex- 
periments here is from the MC Glauber model. 
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FIG. 1: (Color online) Number of charged hadrons at midra- 
pidity in y^sjvjv = 200 GeV Au+Au collisions scaled by the 
number of the participant pairs calculated from the optical 
Glauber model. The data are from the STAR Collaboration 



where T^'' = (e + P)u^'u'' ~ Pg^"" is the energy- 
momentum tensor, is the fluid four-velocity, e is the 
energy density and P is the pressure. As we are interested 
in particle and photon production at midrapidity, we may 
assume that net-baryon density is neghgible. Since the 
particle spectra are approximately flat at midrapidity, we 
can simplify our hydrodynamical equations by assuming 
longitudinal boost-invariance. We use the SHASTA algo- 
rithm [2^, [s^l to solve this (2-|-l)-dimensional numerical 
problem. 

To close the hydrodynamic equations we need an Equa- 
tion of state (EoS), P — P{e). In this paper we study 
the sensitivity of thermal photon production to the EoS, 
by focusing on two different cases. The first case, called 
here " eosQ" , corresponds to the Bag model EoS with 
a first order phase transition [3l|. In eosQ, the high- 
temperature phase with the Bag constant is an ideal 
gas of three flavors of massless quarks and gluons, while 
the low-temperature phase is an ideal gas of all hadronic 
states with m < 2 GeV. These two phases are connected 
with a mixed phase, and the Bag constant is chosen so 
that the critical temperature is T^, = 165 MeV. The sec- 
ond EoS case, which we call "eosL", is adopted from 
Ref. [12]. This EoS is quite similar to the recently con- 
structed lattice EoS " s95-p" |33| , and as discussed in [sl] , 
the hadron spectra and elliptic flow are in practice insen- 
sitive to the differences between eosL and s95-p. 

Since the lattice data suggests that the phase transition 



from the QGP to hadron gas (HG) is not of first order, 
one may consider the eosQ case somewhat unrealistic. 
However, the computation of thermal photon production 
in the phase-transition region requires well-defined QGP 
and hadron-gas fractions, which are available only in the 
eosQ case. With eosL, in the absence of such phase frac- 
tions, there are additional uncertainties in the thermal 
photon calculation related to the QGP and HG emission 
rates. 

Thermal transverse momentum (pt) spectra of 
hadrons are obtained using the Cooper-Frye method [s^l 
where particle emission from a freeze-out hypersurface a 
is calculated with 



dNf 
cPprdy 



f{x,p)p''da^ 



(4) 



where f{x,p) is the momentum distribution function of 
a specific hadron type. We assume the system to de- 
couple at a single constant temperature Tdcc, which is 
fixed so that we get a good agreement with the mea- 
sured pt spectra of pions at RHIC. With eosQ, we have 
Tdoc = 140 MeV, and 160 MeV with eosL . 

After the thermal emission of particles from the freeze- 
out surface is calculated, we take into account the strong 
and electromagnetic 2- and 3-body decays. This treat- 
ment is essential since most of the stable particles in our 
case come from the decays of heavy resonances. 

The Pt spectra of hadrons can be written as a Fourier 
series. 



dN 



1 dN 



d^PTdy TT dpTpdy 



00 

1 + ^ 2t;„ cos(n0)) , (5) 



where cj) is the hadron momentum's azimuthal angle with 
respect to the reaction plane defined by the impact pa- 
rameter. Elliptic flow, V2^ is the second coefficient in this 
series and it can be computed from 



V2{Pt) 



/d0cos(2<^)^^^ 



dN(b) 



(6) 



Correspondingly, the pT-integrated V2 becomes 
/d0cos(20)fg 



V2 



, dN{b) 
d<f)dy 



(7) 



D. Thermal photons emission from the 
hydrodynamical medium 

The Pt spectra of thermal photons can be calculated 
from 



dN-y 

d^PTdy 



d'^xT{E*{x),T{x)), 



(8) 
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FIG. 2: (Color online) The HG-to-QGP ratio of the photon 
emission rates as a function of the photon energy at two differ- 
ent temperatures. Two different HG-rates, R92 [3^ - [40l ] and 
TRG [4JI, are compared. 



case for the R92 rates. 

Elliptic flow for the thermal photons is calculated as 
in Eq. ((5]). Since the thermal photons cannot be distin- 
guished from other direct photons, the elliptic flow from 
thermal photons alone cannot be measured. In what fol- 
lows, we assume that the net contribution to the photonic 
V2 from the other direct photon sources remains small, es- 
pecially since the fragmentation photons with a positive 
V2 should partially cancel the negative f 2 of the photons 
arising from parton-medium interactions [9|]. 

If other components are emitted isotropically we can 
roughly estimate how much they "wash away" the elliptic 
flow coming from thermal photons. The total elliptic flow 
is then 




where dN^^"^ / dp^dy corresponds to the measured pr spec- 
trum of direct photons and ^2'^ is the V2 of thermal pho- 
tons alone. 



III. RESULTS FOR RHIC 
A. Hadron spectra and elliptic flow 



where T{E* ,T) is the Lorentz invariant thermal photon 
emission rate, rf^a; is the volume element and E*{x) — 
p^u^{x) is the photon energy in the fluid's local rest 
frame. For the QGP, we use the emission rate from 
Refs. (ssl-lsSl with Nf = 3 and a running strong coupling 
constant [sl = /3/ ln(8r/A), with /3 = 67r/(33 - 2N}) 
and A = 200 MeV. For the hadron gas, we use two dif- 
ferent emission rates: (i) Those calculated in [3^,[3^ and 
parametrized in [s^, [IqI , which we call " R92" . These 
rates were used in the previously published LHC predic- 
tions [20^] . (ii) The more recent ones from Ref. [4lj which 
account also for the finite size of hadrons through form 
factors. We call these rates "TRG". 

With eosL, which smoothly goes from the QGP to the 
HG phase without specifying their volume fractions, one 
needs to choose how to switch from the QGP to HG 
photon emission rates. For simplicity, we choose to do 
this at a constant temperature Tg but we vary between 
170 and 200 MeV. We label these two cases as "eosL170" 
and "eosL200". 

To illustrate the differences between the R92 and TRG 
emission rates which will be important for the photonic 
V2 results presented in Sec. Ill ahead, we plot in Fig. [2] 
the ratio of the photon emission rates in the HG and QGP 
at two different fixed temperatures. We see that (since 
both HG rates have been divided by the same QGP rate) 
the difference between R92 and TRG is about a factor 
six at large energies. Furthermore, the TRG rates are 
always well below the QGP rates, while this is not the 



First we show the hadronic observables to demonstrate 
that our hydrodynamical description of the bulk QCD- 
medium is reasonable. Figure [3] shows the transverse 
momentum spectra of positively charged pions in differ- 
ent centrality classes in ^/snn — 200 GeV Au+Au colli- 
sions at RHIC. As we can see, we have a good fit to the 
pion spectra below p^ ~ 2 GeV for a very wide range 
of centralities. The integrated elliptic flow of charged 
hadrons from Eq. (O is plotted in Fig. |4] together with 
the data obtained by the STAR Collaboration [43] using 
the 4-particle cumulant and LYZ methods which should 
best reflect the elliptic flow relative to the reaction plane 
defined by the impact parameter. We have a fairly good 
description of the data also here, although the central- 
ity dependence of the computed V2 is not fully repro- 
duced. We expect, however, that fine-tuning the initial 
density profile, invoking the Monte Carlo Glauber model 
and possibly also event-by-event hydrodynamics (see e.g. 
Ref. 143) as well as including viscous effects (see e.g. 
Ref. [45| ) will improve the agreement. These improve- 
ments are, however, beyond the scope of this exploratory 
paper. 

B. Photon spectra 

In order to study how much elliptic flow is washed 
away by other sources of direct photons at RHIC, we 
need to estimate the other components. We do this by 
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FIG. 3: (Color online) The pT spectra of positive pions for 
y/SNN ~ 200 GeV An + Au collisions at RHIC compared 
with the PHENIX data [43 |. The centrality classes are indi- 
cated in the figure and the spectra are scaled by increasing 
powers of 10"^. 
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FIG. 4: (Color online) The integrated V2 for charged hadrons 
in y^sjvjv ~ 200 GeV Au + Au collisions at RHIC compared 
with the STAR data El. 



fitting the measured photon pT spectrum. To study the 
uncertainties due to the chosen fit functions, we use two 
different forms. Our first choice ("fit 1") is an exponential 
combined with a power law function [19| 



/(pt) = Aexp{-pT/T) 



C 



(10) 



where A, T, C, b and n are the fit parameters. This fit 
function is physically motivated by the QCD-like power- 
law behavior at high pT and the thermal-like exponential 
at low Pt- Our alternative choice ("fit 2") is a mere 
power-law function 



f{pT) 



c 



(l+p2,/6)» 



(11) 



In these fits, we use the photon data from PHENIX 
Collaboration. The older data sets [13, [3 have large 
error bars at low px but in the more recent low-p^ data 
the error bars are much smaller. Unfortunately, the 
centrality classes in these measurements differ from each 
other. For our fits shown in Fig. [SI we have simply com- 
bined the 0-20% (20-40%) centrality data from Ref. Jll 
with the 10-20% (30-40%) centrality data from Ref. [l8j . 
In our fits we have included all datapoints from the above 
sets. 

With the fit 1, we first find the parameters b and n by 
fitting the measured photon px spectra in p-|-p-collisions 
using the PHENIX data [o, 111 shown in Fig. E^a). 
Then for the Au-|-Au case, keeping the high-pT slope- 
parameters b and n fixed, we find A, T and C by fitting 
the PHENIX data [H, for the two centraHty classes 
shown in Fig. [5j For the fit 2, we use the same data sets. 
The best fit parameters obtained for the power law fits 
are listed in Table HIl and the parameters for the fit 1 can 
be found in Table Hill The fits 1 and 2 have equally small 



values at both centralities. 



Power law fit 





C [GeV""] 


b [GeV^] 


n 


p+p 


3.29- 


10-^ 


4.37 


10-^ 


3.09 


Au-hAu 0-20% 


2.16 


10" 


5.24 


10"^ 


3.35 


Au-I-Au 20-40% 


9.33 


10^^ 


6.35 


lo-'* 


3.52 



TABLE II: The parameters obtained for the power law fits 2. 



Exponential -I- power law fit 





A [GeV 


] T [GeV] 


C 


Au-I-Au 0-20% 


85.4 


0.212 


4.96 


Au-hAu 20-40% 


30.7 


0.218 


1.18 



TABLE HI: The parameters for the exponential 4- power law 
fits 1. In this case, n and 6 are obtained from Table HIl 

In Fig. EJb) we have replotted the low- and mid-pT 
region from Fig. Ej^a), and shown our thermal photon 
results obtained with eosL (eosQ) using the TRG (R92) 
rates in the HG phase. For clarity, we have plotted the 
eosL results only for Ts = 170 MeV. If we do the switch of 
the emission rate at T, = 200 MeV, we get 30% (10%) less 
photons at pT — 1(2) GeV, because the TRG emission 
rate is smaller than the QGP emission rate, as shown in 
Fig. [2 

Our thermal photon results shown in Fig.[5Jb) differ by 
a factor of two at high px- Some of this difference comes 
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FIG. 5; (Color online) The fits 1 and 2 to the measured photon spectra for two diflferent centralities in -y/ijvjv = 200 GeV 
Au+Au collisions and in p+p collisions at ^/s — 200 GeV. The 20-40% centrality-class spectra are divided by 100 and the p+p 
case by 1000. The data are from the PHENIX Collaboration [l7l - [l^ . Our thermal photon results are shown in the panel (b). 



from the small difference in the initial temperature pro- 
files which in our case are obtained from the fixed initial 
entropy density through the EoS. However, a more domi- 
nant effect is the different mapping of the energy density 
to the temperature in eosQ and eosL. The difference in 
the actual temperature in the two cases is not large but 
the exponential temperature-dependence in the emission 
rates magnifies the effect considerably. 

From the previous photon studies, see e.g. Refs. [Ill, 
[T^ . we know that photons from HG are contributing 
mostly at small px- In Fig. [5] the difference between the 
eosQ+R92 and eosL-l-TRG results shrinks down at low 
Pt since the R92 HG emission rate is larger than that 
in the TRG rates and since with eosQ the HG volume 
becomes larger than with eosL. When we use the TRG 
rates and eosL, only 3% of the photons come from the 
HG at Pt = 1 GeV. With eosQ and the R92 rates about 
50% of the photons originate from HG at the same px- 

We also note that in the \ow-pT region our ther- 
mal photon results very clearly undershoot the latest 
PHENIX data. We have checked that changing the 
freeze-out temperature to Tdoc = 120 MeV gives only 
a negligible improvement. This feature is typical to al- 
most all hydrodynamical calculations as can been seen 
e.g. from Fig. 43 in Ref. [46*1 . 

At ~ 3 GeV the obtained thermal photon emission 
is almost enough to match the fit 1 at both centralities 
if eosL is considered. This suggests that we may have a 
window for thermal photon dominance at this pr- How- 
ever, if we compare with the fit 2 there is always at least 
a factor of two difference. Event-by-event fluctuations in 



the initial state, however, have been shown to increase 
the thermal emission at pt > 2 GeV [47|, and thus we 
should indeed have a better chance to have a region where 
the direct photon px spectrum, and consequently also the 
photon V2 , at RHIC is entirely dominated by the thermal 
emission. 

C. Photon elliptic flow 

In Fig. inja) we have plotted the elliptic flow of the 
thermal photons using both eosQ and cosL with the R92 
rates. Panel b shows the same calculations but with the 
TRG rates. Unlike for hadrons, the thermal photon el- 
liptic flow starts to decrease quickly above pT ^ 2 GeV. 
The reason for this is that practically all high pT photons 
are emitted nearly isotropically in the beginning of the 
evolution (see e.g. Fig. 3 in Ref. [13]), when the hydro- 
dynamical flow effects are very small. Since the photon 
emission is dominated by the early times the thermal 
photon elliptic flow is clearly smaller than the hadronic 
W2(pt), which probes the flow anisotropy only on the 
freeze-out surface. 

From Fig.[6ja) we see that the larger switching temper- 
ature Ts in eosL only moves the elliptic-flow peak towards 
higher pr- However, as shown in Fig.lS^b), with the TRG 
rates there is factor of two difference in the maximum 
value between the eosL170 and eosL200 cases. This sys- 
tematics can be deduced from Fig. [5] and hydrodynam- 
ical evolution as follows: When the system reaches the 
cross-over region near the QCD phase transition, there 
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FIG. 6: (Color online) Elliptic flow of thermal photons in Au+Au collisions at ^/snn = 200 GeV. 



has already been a significant anisotropy developed in the 
transverse flow which is directly reflected to the thermal 
photon V2- Thus, if the photon emission is increased (de- 
creased) in this region, the thermal photon V2 increases 
(decreases). As seen in Fig. [21 with the TRG rates the 
QGP emission rate is larger than the HG rate. Hence 
increasing the switching temperature decreases the emis- 
sion and thus the V2- This effect is seen in Fig. [6jb). 
With the R92 rates in Fig. El^a) the situation is slightly 
different, as in the cross-over region T = 170. ..200 MeV 
the QGP emission rate is larger than the HG rate at small 
energies and vice versa at high energies. Thus, at \ow-pT 
in Fig. inja) the situation is similar to panel b (i.e. V2 is 
larger for eosL170 than for eosL200). At > 2 GeV, 
the increase of the switching temperature now increases 
the total emission and thus making V2 larger for eosL200 
than for eosL170. 

We also notice from the eosQ results in Fig. |6] that the 
maximum V2 decreases by a factor 2 when replacing the 
R92 rates by the TRG rates. Since the QGP rates in both 
cases are the same, this signals to us that the hadron gas 
indeed plays an important role in generating the thermal 
photon V2 in the eosQ case. To quantify this statement, 
we have plotted in Fig. [7] the fraction of photon V2{pt) 
coming from the HG phase. We define v^'^ as 



,.HG 



J dip cos(20) 



dN" 



dp^d(pdy 



I' 



dp'i^d(f)dy 



(12) 



i.e. relative to to all thermal photons. We see that the 
photon V2 can be mostly from the HG (eosQ with R92) or 
mostly from the QGP (eosL170 with TRG), or between 
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FIG. 7: (Color online) The contribution to the thermal pho- 
ton elliptic flow from the hadron gas in Au+Au collisions at 
^/sJTn = 200 GeV. 



these extremes (eosQ-|-TRG and eosL170-|-R92). Thus, 
both the EoS and the HG emission rate have a big effect 
on where the thermal photon V2 originates from. 

Figure [51 illustrates how much elliptic fiow of thermal 
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photons is washed away if we include other direct photon 
components assuming that they are produced isotropi- 
caUy. We can see that the final photon V2 , obtained from 
Eq. ^ based on the fits 1 and 2, is clearly smaller than 
the thermal one, and also that the different fit functions 
modify the place and shape of the peak, keeping however 
the maximum V2 roughly the same. 
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FIG. 8: (Color online) Elliptic flow of thermal and direct 
photons in Au+Au collisions at y^ijvjv = 200 GeV. 

Finally, in Fig. [HI we have plotted both the thermal 
photon and the full direct photon elliptic fiow in 0-20% 
and 20-40% centrality classes for the eosQ-|-R92 and 
eosL170-|-TRG cases. The latter can be considered as 
a state of the art calculation in that a realistic EoS and 
latest rates are utilized. In the eosQ-|-R92 case, the el- 
liptic flow is as large as it can be in our approach. For 
both cases, the fit 1 is used to estimate how the other 
components reduce the elliptic fiow. As seen in Fig. [5l 
the thermal photon yield is smaller in the eosQ case, and 
hence the other direct photon components wash away 
more of the elliptic flow in the eosQ case than in eosL170 
case. 



IV. PREDICTION FOR THE LHC 

Next, we extrapolate our hydrodynamical modeling to 
the .JJWn = 2.76 TeV Pb+Pb colhsions at the LHC. 
We choose the same sWN initial density profile and de- 
coupling temperatures as at RHIC, and, as explained 
in Sec. II. B, use the measured charged-hadron multi- 
plicity [25} to fix the initial entropy and initial time 
through the EKRT-model (for details, see Ref. [1^). As 
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FIG. 9: (Color online) Elliptic flow of direct photons in 
Au+Au collisions at y^ijvjv ~ 200 GeV. 



seen in Fig. |TU1 a reasonable agreement with the mea- 
sured charged-hadron pT spectrum (4^ follows up to 
Pt ~ 4 GeV. This ensures that our thermal photon cal- 
culations are meaningful also at the LHC energies. 

In Fig. |TT]we have plotted our prediction for the ther- 
mal photon Pt spectrum in the ^sjvAf = 2.76 TeV 
Pb+Pb collisions at the LHC. The bands shown are de- 
fined by the cases eosL170+R92 and eosQ+TRG, which 
give the largest and smallest yields, correspondingly. 
We see that the uncertainty coming from the EoS and 
from the HG emission rates is at largest of the order 
of 40%. We note that these predictions are qualita- 
tively quite similar to the predictions given in [20l [2]| 
for y^IWJJ = 5.5 TeV. 



9 




FIG. 10: (Color online) Transverse momentum spectra of FIG. 12: (Color online) Elliptic flow of thermal photons in 
charged hadrons in 0-5 % most central Pb+Pb collisions at Pb+Pb collisions at ^/snn = 2.76 TeV. 
^si^ = 2.76 TeV. Data from ALICE Collaboration 




FIG. 13: (Color online) The contribution to the thermal pho- 
FIG. 11: (Color online) Transverse momentum spectra of ton elliptic flow from the hadron gas in Pb+Pb collisions at 
thermal photons in Pb+Pb collisions at ^sjvjv ~ 2.76 TeV. ^/sNN = 2.76 TeV. 



Since currently we do not have the measured total di- 
rect photon spectrum at the LHC available yet (like we 
had at RHIC), we can consider only the thermal pho- 
ton elliptic flow here. It is, however, very interesting to 



compare the thermal photon elliptic flow with the RHIC 
results. We have plotted the obtained thermal photon V2 
in Fig. [12] for 20-40% central collisions. We can see that, 
similarly to the hadronic case [i^, the thermal photon 
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elliptic flow is very similar at RHIC and LHC. This is a 
non-trivial result since the temperature-range, flow-range 
as well as the volume factors for photon emission (from 
the QGP in particular) are larger at the LHC than at 
RHIC, and, as discussed in ^] also the flow asymmetry 
near the phase transition region is larger at the LHC. 
Then, when going from RHIC to LHC, in order to arrive 
at a similar V2 in both cases, the increased flow asym- 
metry in the numerator of Eq. |6] is compensated by the 
increased photon yields in the denominator of Eq. [5] In 
Fig. [13] we again plot the hadronic fraction of V2{pt) for 
the same cases as in Fig. [T] The flgure shows that vf^ is 
again very close to the corresponding fraction at RHIC. 

V. DISCUSSION 

We have considered the sensitivity of thermal photon 
production to the EoS and emission rates in heavy-ion 
collisions at RHIC and LHC. We have compared the ob- 
tained thermal photon yields with the PHENIX measure- 
ments, and shown that in the window 2 < ^ 3 GeV 
the thermal contribution, computed with a realistic EoS 
(eosL) and latest emission rates (TRG), is reasonably 
close to the data. Like in most previous hydrodynami- 
cal studies, in the region pT ^ 1 GeV, however, we get a 
clearly smaller yield than what is measured most recently. 
We have shown that around ^ 2 GeV the thermal 
photon elliptic flow peaks at a fairly large value, 5 % 
in the 20-40% centrality class with eosL-|-TRG, but also 
that the possible other components may wash even half 
of this away. We emphasize, however, that the amount of 
V2 wash-out depends on the thermal photon contribution 
relative to the other components. Thermal photon pro- 
duction near px ^ 2 GeV at RHIC can be expected to in- 



crease further once the event-by-event QCD-matter den- 
sity fluctuations are accounted for [J^l , in which case the 
thermal photon production can become dominant and 
the V2 wash-out to decrease or even vanish. 

Constraining our hydrodynamical modeling with the 
measured charged-hadron spectrum in ^/s'nn — 2.76 TeV 
Pb-t-Pb collisions at the LHC, we have predicted the ther- 
mal photon pt spectra and V2- According to our results, 
elliptic flow of thermal photons at the LHC and RHIC are 
very similar in the few-px region. For the determination 
of a possible thermal photon window, and consequently 
thermal photon V2, it will be extremely interesting to see 
the direct photon data at the LHC. 

Next, one should consider the eff'ects of event- by- event 
density fluctuations [3, [50l - [&3 | on thermal photon el- 
liptic flow. On the theoretical side, one would need a 
better understanding of how the degrees of freedom in 
the QGP would be best accounted for when computing 
thermal photon production, as well as a better control 
over the photon emission in the phase transition region. 
Also the dissipative hydrodynamical effects to thermal 
photon production should be studied further, so far only 
the very first steps into this direction have been taken, 
see Ref. [s^. 
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